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Abstract 

We examine a strategy for using neutral current measurements in long- 
baseline neutrino oscillation experiments to put limits on the existence of more 
than three light, active neutrinos. We determine the relative contributions of 
statistics, cross section uncertainties, event misidentification and other system- 
atic errors to the overall uncertainty of these measurements. As specific case 
studies, we make simulations of beams and detectors that are like the K2K, 
T2K, and MINOS experiments. We find that the neutral current cross section 
uncertainty and contamination of the neutral current signal by charge current 
events allow a sensitivity for determining the presence of sterile neutinos at the 
0.10-0.15 level in probablility. 
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1 Introduction 



In recent years a series of exciting experimental results have shown that neutrinos have 
finite masses and mixings. For a recent review of the status see Ref . . Solar neutrino 
and atmospheric neutrino results indicate that all three known neutrino flavors (e, 
fx, t) participate in neutrino mixing, and hence neutrino oscillations. Consequently, 
the standard framework to describe the experimental results and analyse neutrino 
oscillation data is that of three-flavor mixing in which the three flavor eigenstates are 
related to three mass eigenstates by a 3 x 3 mixing matrix[2|. The positive signal 
for P M — > v e oscillations from the LSND experiment [2| challenges the three-flavor 
mixing paradigm 4J. However, the neutrino oscillation interpretation of the LSND 
observations is yet to be confirmed. Independent of whether or not the LSND results 
are confirmed by MiniBooNE[5], the three-flavor mixing framework deserves further 
experimental scrutiny in the coming years. Much of the focus on future experiments so 
far has been directed to the determination of the 3-mixing angles and the CP-violating 
phase with long-baseline or oscillation experiments .6] and reactor experiments jZj. Of 
interest in this paper is the measurement of the neutral current, which could allow 
tests of the unitarity of the 3x3 mixing matrix and thus indirectly probe the existence 
of sterile neutrinos. 

In a three-flavor neutrino model, the sum of the oscillation probabilities 
Y2 y =e ii t P '( u x ~^ v y) is unity. If there are more than three light neutrinos, we know 
from measurements of the invisible width of the Z jH] that the additional neutrinos 
must be sterile. If additional light neutrinos mix with the three known flavors we can 
expect a non-zero oscillation probability to sterile neutrinos, P[y x — > v s ) ^ 0. To test 
the three-flavor neutrino-mixing paradigm it is important to search for a sterile neu- 
trino component within the neutrino flux from natural and manmade sources. Since 
sterile neutrinos have no strong or electroweak interactions, they cannot be detected 
directly. However, neutral current (NC) measurements allow J2 y=e u t P{ v x. —* v y) 
to be determined which, by probability conservation, is equal to 1 — P[y x — > u s ). 
Therefore, in principle a NC measurement alone is sufficient to determine P{u x — > u s ). 
However, in a realistic detector misidentifications of CC and NC events, together with 
systematic uncertainties on the relevent neutrino interaction cross sections, compli- 
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cate the analysis. 

In this paper we study the use of NC measurements to determine limits on the 
sterile neutrino content in long-baseline neutrino oscillation experiments. First we 
consider the sensitivity to the sterile content that might be obtained in a K2K-like 9J, 
T2K-like[m], and MINOS-likejll] experiment with a "perfect" detector and "perfect" 
beam if there are no systematic uncertainties. We then consider the impact on the 
sensitivity of event misidentification and systematic uncertainties. Our study is based 
on a simple simulation of the long-baseline neutrino beams, neutrino interactions |12j . 
and detector responses. We present our results versus event rates and the size of the 
cross section uncertainty in order to show the dependence on these quantities. 

2 Using NC data to determine sterile content 
2.1 Formalism 

The present and proposed long-baseline neutrino oscillation experiments exploit con- 
ventional neutrino beams that are produced by the decays of charged pions in a long 
channel. This produces a beam which is initially almost entirely v^. Kaons and 
muons decaying in the channel introduce a small (typically ~ 1%) v e component in 
the neutrino beam. As the neutrino beam travels towards a distant detector its flavor 
content will evolve. In our analysis we will consider three active neutrinos (i/ e , u^, 
v T ) and one sterile neutrino (z/ s ), with oscillation probabilities P[y^ — > v x ) = P/xx, 
x = e, fi, t, s. We begin by considering an oscillation experiment that has an initially 
pure beam with well known neutrino spectrum and flux, and a detector with per- 
fect identification of the produced events. Then the event rates at the far detector 
will be 






(2) 
(3) 
(4) 
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where N® is the predicted number of CC interactions in the detector in the absence 
of oscillations and the a x denote the interaction cross sections (a e , <r M , a T ) for (u e , 
v T ) CC interactions and ctnc for v x NC interactions. 

In an ideal experiment A^c determines P MS and Acc determines P^. In practice 
the presence of v e and v T CC events complicates the analysis if these events are not 
distinguished from NC events. In that circumstance the NC events provide a measure 
of 1 — P MS + e e P Me + e T P^ T , where the factors e e and e r reflect the contaminations. 

Probability conservation (P Me + P w + P MT = 1) can be used to eliminate P MT or 
P Me , but not both. If the beam energy is below the threshold for r production or the 
probability P Me is small and can be neglected, then P^ can still be determined. How- 
ever, for a realistic detector with particle misidentifications and/or a u e component 
at the far detector that cannot be neglected, the problem of determining P^ s can be 
complex but still solvable, as we shall discuss. 

In general, let the probability that an event of type x (NC, v e CC, CC, or 
v T CC) be identified in the detector as an event of type y be given by ( xy , where 
x, y = (NC, e, fx, r) (note that ( xx is the efficiency for detecting an event of type 
x). If N® is the predicted number of CC interactions in the detector in the ab- 
sence of oscillations, then after including oscillations, detector efficiencies and mis- 
identifications, and integrating over the energy dependence, the number of measured 
events of type y will be: 



where the interaction cross sections for v e CC, CC, v T CC, and NC events are 
given by a e , a^, o T and ctnc, respectively. 

2.2 Ignore z/ e 's 

We consider first the situation in which the v e component in the beam at the far 
detector is so small that v e CC interactions can be neglected. In this case we let P^ e — > 
(it is known to be small, at the 5% level or less from the CHOOZ experiment |13j). 



N, 



y 
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It is convenient to define the following two ratios, 



N. 



RnC = 7 N< atO ~ (1 ~~ + fvhNcPfjtn + f^NcPfir , (6) 

(.NCNC 1 ^NC 



3> 



where 



and 



iV° = o-,JVj/a„ , (8) 



— CxyCx/ Cyy^y (9) 

is a normalized misidentification factor that gives the fraction of events identified as 
being of type y that are really events of type x. Measuring R NC and is sufficient 
for deducing P^ s (and P^)- The analysis depends on whether or not we are above 
the v T CC interaction threshold, i.e., whether or not there are v T CC events produced 
in the detector. 



2.2.1 Below r threshold 

For neutrino energies below the r threshold o T = and f T j = 0. In this case we can 
invert Eqs. El and [7| to obtain 

MA* ~~ 1 f r ' 
J- — J^NCJNC^ 

p _ i R NC ~ Rpfn,NO ny) 

A _ JH,NCJNC,iJ, 

Adding uncertainties in quadrature we get 



SPw = 1 ; — ; ; > ( 12 ) 





+ fNcA SR ™) 2 


i - 


fn,NcfNC,iJ, 


\J (SRnc 


) 2 + KncW 



SP,s = 1 : ; - f , (13) 

1 — Jn,NCjNC,fi 



where in the limit of Gaussian statistical uncertainties 



5R 3 =R JX j^ + e], (14) 



and 



(15) 



The first term in each <5i?j is the usual statistical uncertainty, the second comes from 
the normalization uncertainty (flux and cross section). 

Note that the normalized mis-identification factors f^c and Jnch will be sensi- 
tive to the neutrino energy spectrum and the detector technology, and therefore must 
be evaluated for each experimental setup. Most of the mis-identification terms are 
suppressed by f 2 ; if / < 0.1 then f 2 < 0.01. If the experimental setup is such that 
we can ignore all terms of order f 2 , (see Table 1) we have 

1 — RnC + f/j,,NcRp 



(16) 



SP^s SR NC 
which measures the significance of the deviation of P^ s from zero. 

For a perfect detector that can identify each event correctly, f x>y = 5 xy . In this 
limit P^ s = 1 — Rnc and 



P P 

1 us 1 us 



;i7) 



This ratio depends only on P^ s , the experimental statistics, and the systematic un- 
certainty on the NC measurement. Thus, Eq. El defines the maximum sensitivity 
that is in principle achievable for a given N^ c and e^c- 

2.2.2 Above r threshold 

If the neutrino energy is above the r threshold and there is not a clean signature for v T 
CC events, we can still deduce P w and P MS by using the identity P w + P MT + P^ s = 1 
to eliminate P MT in Eqs. Eland[7| (we are still assuming P^ e = 0), which gives 

p _ P M ( 1 + fr,Nc) + RncUnC^ + /t,/J 

1 + fr,NC — fr,n + fr,NcfNC,n ~ fn,Nc{fNC,/j. + fr,p) ' 

p _ j Rnc(1 - fr,n) - RjJJ^NC - fr,Nc) ^gs 

1 + fr,NC — fr,n + fr,NcfNC,n ~ fp,Nc{fNC,iJ, + /r,/i) 

If no other process contaminates the CC events (i.e., fj tfl = as appears to be the 
case for a MINOS-like experiment; see Sec. EJ), then 

P/J.S ^ 1 + fr,NC ~ RNC + Rftif^NC ~ fr,Nc) 



SP»s ^(6R NC ) 2 + (U,nc - fr,Nc) 2 (5R,) 2 
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(20) 



For a perfect detector, P MS /5P MS is again given by Eq. [T7| 



2.3 Do not ignore z/ e 's 

If the v e CC interaction rate in the far detector is not negligible (which could be the 
case if sin 2 2# 13 is near its upper bound and we want to push the uncertainty in the 
measurement of down to the few per cent level), then we need three measurements 
to be able to solve for all of the probabilities. The potential measurables are 

N NC 

RNC = 7 TTq - = (1 ~~ Pfis) + fn,NcP^n + fe,NcPfie + fr^cPfir , (21) 

^NCNC^NC 

^ ^ 7 NO = f NC 'V^ ~ ^A») fe,/J,Pfie + fr^Pfir > (22) 

W iV At 

CeeN° e 

and, if we are above the v T CC interaction threshold, 

Rt = ^ = /iVC,<r(l — Pfis) + f^rP^ + fe,rPfie + P/jt ■ (24) 



R £ = 7 — iTrn ~~ /iVC,e(l — Pfis) + ffi^Pfifj, + Pfie + fr^Pfir , (23) 



2.3.1 Below r threshold 

Below the v T CC threshold energy the three measurements must be P M , Rnc, and 
R e . Then / T j = 0, the P^ T terms drop out, and we can invert Eqs. l2"Tll2*3l to obtain 



p R[i(l — fe,NcfNC,e) ~ RNcjfNC,^ ~ JiVfie Je,/J ~ Re(fe,fi ~ fe,NcfNC,^) / of -x 

W — 1 — / ' 

P e (l — f^NcfNCtt) — RNc{fNC,e ~ fNC,^ffj,,e) ~ Rfi(f^,e ~ j 'u,NC j ? NC,e) /0 „ x 
P »e = YZTf ' ( 26 ) 

p _ , PjVCjl ~ fft,efe,n) — Rft(ffi,NC ~ fjj,,efe,Nc) ~ Re{fe,NC ~ fe,fif[i,Nc) (r>7 \ 
*fJ.s — i I — f ' 

Where / = f^NcfNC^ + fe,NcfNC,e + ffi,efe,[i ~ f^,NcfNC,efe,^ ~ f NC,fif fi,ef e,NC ■ The 

calculation of the £P's is straightforward; each P term has a statistical and systematic 
uncertainty given by Eq. El 

Note that for an idealized detector in which no other processes significantly con- 
taminate v e CC events (i.e., / J)e ~ 0) and v e CC events do not contaminate CC 
events (i.e., / 6iAt ~ 0), then P Me = P e . Since P Me is small (of order 0.1 or less, as 
indicated by current oscillation limits), eliminating terms of order f 2 and fP^e in this 
case will recover the situation where we ignored v e (i.e., Eq. H7)|) . 
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2.3.2 Above r threshold, no r measurement 

For energies above the v T CC interaction threshold the P^ T terms do not drop out of 
Eqs. If we do not have the means to measure v T CC events but can measure 

v e CC events, then we can use probability conservation to eliminate P Mr , giving 

RNC = (1 — PfiajO- + fr,Nc) + P^{ffi,NC ~ fr,Nc) + Pfie{fe,NC ~ fr,Nc) 5 (28) 
R, = (l-P,s)UNC,» + U)+P^(l-U)+P»e(fe,»-U), (29) 
Re = (l-^fee + W+^^-W+^fl-g. (30) 

The general solution for the probabilities is somewhat messy, but if we assume that 
no other processes contaminate the CC signal (i.e., fj 7fl ~ 0) and the v e CC events 
do not contaminate the other signals (/ ej - ~ 0), (see Sec. |3J), then P m = R^ and we 
can invert Eqs. EHl and EH to obtain 

p _ ^ _ Rnc(1 ~ fr.e) + R e fr,NC + Ry [fr,Nc(^ ~ f^e) ~ f^,Nc(^ ~ /r,e)] /g-Q 

1 + /r,JVc(l + /iVC.e) — fr,e 

The calculation of SP^g is straightforward. 

If no other processes contaminate the v e CC signal (i.e., fj e — > 0), then P^ e = R e 
and we obtain 

Pfis 1 + fr.NC ~ RNC + Rn(fn,NC ~ fr.Nc) ~ Refr,NC 



6P » S J(5R NC ) 2 + (U,NC - fr,Nc) 2 (5R,) 2 + f* NC ($Re) 2 



(32) 



2.3.3 Above r threshold with a r measurement 



If R T is also measured, in addition to R e , then there are four measurements (P M , Rnc, 
R e and R T ), but there are only three independent quantities (since P w + P Me + P MT + 
P^s — 1)- One possible approach would be to assume that P^ s is independent of the 
other probabilities and use these four measurements to test probability conservation. 
We do not pursue this option here. Instead, we use probability conservation to elimi- 
nate one of the probabilities and use three of the four measurements to determine P MS 
(the fourth measurement could be used to check probability conservation afterwards). 
Since Piv^ — > u T ) is most likely much larger than Piv^ —>■ v e ) in the L/E regime 
we are considering, we use R T as the third measurement (along with P M and Rnc)- 
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Then the appropriate formulas for the measurables Rncj R^ an d Rr can be found by 
the interchange r e in Eqs. 

3 Detector simulations 

We wish to explore how well in principle a neutrino three-flavor unitarity test can be 
performed with a given muon-neutrino beam as a function of dataset size, and study 
which systematic uncertainties are likely to be important, and their impact. 

We consider first a "perfect" experiment in which the sensitivity of the unitarity 
test is determined only by the statistical uncertainties, calculated using a parameter- 
ization of the known beam flux and spectrum, together with a simulation of neutrino 
interactions in the detector. An event simulation is used to determine the relevent de- 
tection efficiencies and misidentification factors. We use the NEUGEN Monte Carlo 
code ^2] to simulate neutrino interactions in the detector. Events are classified as 
Up CC, u e CC, or NC. In practice the requirements used to identify events of a given 
type will depend upon the detector technology. For example, for a water cherenkov 
detector in our simple analysis we will define a u e CC event candidate as an event 
with an electron candidate above threshold. An electron candidate is either a real 
electron or a tt° with an energy exceeding 1 GeV (in which case the two daughter 
photons from the high energy 7T° produce cherenkov rings that overlap in the detector 
and cannot be distinguished from a single electromagnetically showering particle). 
A NC event candidate would be an event containing a tt° candidate but no muon 
candidate, where a tt° candidate has two e-like rings above threshold (which come 
from a tt° with energy less than 1 GeV). The definition of CC and NC events can of 
course be varied, and then tuned to give favorable values for the signal efficiencies 
and mis-identification factors. Examples are shown in Table 1. 

3.1 K2K-like and T2K-like Experiments 

To identify the most important systematic uncertainties it is useful to compare the 
sensitivity of our "perfect experiment" with that of a realistic experiment. We begin 
with the K2K experiment. K2K uses a beam from the KEK laboratory in Japan. The 
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neutrinos in the KEK beam have a mean energy of 1.3 GeV [9 3 , and the neutrinos 
travel 250 km to the Super-K water Cerenkov detector. A new experiment T2K is 
being planned that will exploit a more intense neutrino source that is presently under 
construction at Tokai, Japan. T2K will also use the Super-K detector, but with 
a slightly longer baseline (300 km) and narrow-band beam with an axis displaced 
slightly from pointing directly at the far detector (an "off-axis" beam). The real 
experimental sensitivities of the K2K and T2K experiments can only be determined 
by the experimental collaborations. In the following we use the NEUGEN Monte 
Carlo program to simulate neutrino interactions together with a simple model for 
the response of a Super-K-like detector. Although this is inadequate to precisely 
predict the real K2K and T2K sensitivities, it does enable us to identify the dominant 
sources of systematic uncertainties, and hence explore how the experimental results 
will depend upon the sizes of these systematics. We use the following parameterization 
of a Super-K-like detector response: 

(a) A threshold of 197 MeV/c for the detection and measurement of muons |14j . 
and 100 MeV/c for electrons and 7r°'s. These thresholds approximate those 
used for the atmospheric neutrino analysis of Super-K E] . 

(b) Energy resolutions given by [TH] 

AE rms 0.025 , , 

— ^£ = 0.005 + ^^, 33 
E \/{E) 

for electrons and 7r°'s and 

^^ = 0.03, (34) 
V 

for charged pions and muons. 

In addition, we use a parametrization of the spectra for the K2K and T2K neutrino 
beams. 

In our analysis we will use only simulated events with visible energy greater than 
0.1 GeV. For our "basic" signals we define a CC event candidate as an event with a 
single muon-like ring, a v e CC event candidate as an event with a single e-like ring, and 
a NC event candidate as an event with two e-like rings, which are assumed to be two 
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Table 1: Signal efficiencies (Cjj) an d normalized mis-identification factors (fij) in 
selected long-baseline experiments. 



Experiment 


j (channel) 


Signal 


Cjj 


JNC,j 


f 


f 

JG,j 


Jr,j 


rvzlv-liKC 


INC 


two e-like, no /t-like 


U.o91 




U.Uuo 


U.Uoz 




(basic) 




one /i-like, no e-like 


U.DZU 


U.Uo f 




U.UUU i 






e 


one e-like, no /i-like 


n a 07 
u.4y / 


U.UUo 


U.UUU4 






rvzlv-llKe 




even e-like, no /i-like 


U.4o / 




U.U < o 


U.UoU 




(aggressive) 


/i 


even e-like, > 1 /t-like 


u.yoy 


U.UoO 




U.UUo 






e 


odd e-like 


u.yyo 


U.UUZ 


U.UUo 






rvZlV-llKe 


l\C 


even e-like, no /i-like 


fl AClA 

U.4y4 




U.Uol 


n m 1 
U.U11 




[ Vjdubhlclll 


M 


even e-like, ^ 1 /i-like 


u.yy^ 


U.U / O 




U.UUU / 




b6<im) 


6 


oqci e-iiKe 


n QQQ 


u.uuuo 








T2K-like 


NC 


even e-like, no /i-like 


0.420 




0.25 


0.006 




(Gaussian 


M 


even e-like, > 1 /t-like 


0.988 


0.036 




0.0014 




beam) 


e 


odd e-like 


0.944 


0.00002 


0.00001 






MINOS-like 


NC 


no (i > 1 GeV 


1.000 




0.903 




0.429 


P^ = 


M 


any /t > 1 GeV 


0.749 












MINOS-like 


NC 


no e, /i, or 7 


0.520 




1.067 


0.005 


0.347 




M 


any /t > 1 GeV 


0.749 















e 


no /i > 1 GeV, > 1 e or 7 


0.999 


0.125 


0.090 




0.064 



photons from a single it decay. Given these definitions, the detector efficiencies and 
mis-identification factors determined from our simulations are listed in Table ^ As 
shown in the table, the efficiencies Qj are of order one-half, and there is no significant 
contamination of one signal by another due to mis-identification. Also shown are the 
results of a more aggressive signal definition, where a simulated event with an odd 
number of e-like rings is labeled as a v e CC event candidate, and the remaining events 
(those with an even number of e-like rings) are labeled as CC event candidates 
if they have one or more /i-like rings or NC if they do not. In this more aggressive 
scenario no events are discarded, i.e., all events were used for one of the targeted 
signals. Although some of the misidentification factors are slightly larger for the 
aggressive scenario, overall they are not greatly changed, while there is a significant 
improvement in the efficiencies for the CC events. 
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To investigate whether our analysis is sensitive to the assumed details of the neu- 
trino spectrum we have repeated the calculation of efficiencies and misidentification 
factors for a K2K-like experiment with a beam that has the same average energy and 
beam spread as the K2K beam, but with a Gaussian energy spectrum (no long high- 
energy tail). For the Gaussian beam, the misidentification factors involving v e were 
greatly reduced (since backgrounds from the high-energy tail are now suppressed), 
but f^NC and /Wc> were only slightly affected. Since f^NC is the dominant / factor 
for a K2K-like experiment, we conclude that our results are not very sensitive to the 
detailed beam spectrum we assume. 

We now consider a T2K-like experiment, where we have used a beam spectrum 
that corresponds to a detector 2 degrees off-axis. The resulting mis-identification 
factors for a T2K-like experiment are shown in Table ^ All of the misidentification 
factors are reduced except for f^NC, which is now 0.25. Therefore, in both the K2K- 
like and T2K-like experiments, the most important contamination is CC events 
being mis-identified as NC events. 

3.2 A MINOS-like experiment 

The MINOS experiment is a long-baseline oscillation exeriment that will use a neu- 
trino beam from the Fermilab Main Injector and an iron-scintillator sampling calorime- 
ter 730 km away in Minnesota. MINOS is expected to begin data taking early in 2005. 
With a beam energy that is about a factor of three higher than the KEK beam, and 
a detector that is very different from the water cerenkov detector used by K2K and 
T2K, the efficiencies and misidentification factors for MINOS will be very different 
than those for the experiments in Japan. To compute the numbers given in Table 
we have once again used a parametrization of the neutrino beam spectrum, the NEU- 
GEN Monte Carlo Program to simulate neutrino interactions in an iron detector, and 
a simple parametrization of the response of a MINOS-like detector. In particular we 
assume: 

(a) An energy threshold of 50 MeV for the detection and measurement of electrons, 
and charged and neutral pions, and a threshold of 1 GeV for the identification 
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and measurement of muons. Note that the MINOS detector is expected to 
be able to determine the charge and measure the momenta of muons from 
0.5 GeV/c to 100 GeV/c, and to distinguish CC events from NC events if 
the muons have momenta exceeding about 1 GeV/c [To] . 

(b) Energy resolutions given by 



for electrons and 7r°'s, 



for charged pions, and 



for muons. Note that in practice the muon energy resolution for the MINOS 
experiment is expected to be somewhat better (worse) than described by Eq. |S7| 
if the muon ranges out (does not range out) in the detector. We found that 
Aprms/p values as high as 0.10 do not appreciably change our results. 

As shown in the table, for a MINOS-like experiment there is a very large contam- 
ination of the NC channel by CC events, and mis-identification of v T CC events 
as NC events is also significant. The efficiency for identifying NC events is about 
one-half, similar to the K2K-like and T2K-like experiments 1 . 

1 Although to first order for MINOS all events with an electron or photon candidate will be 
classified as NC events, there are three independent probabilities, and it is necessary to extract a 
separate v e signal, in addition to and NC signals, to be able to solve for all of the probabilities. 
Hence we must try to select genuine v e interactions from the large NC background. In the table we 
also show mis-identification factors when all non-/Lt events are classified as NC events, which could 
be used when P^ e is very small. 



/S.E rms 0.23 /oi-n 

(35) 



E 



\f{E) ' 



/\E rms 0.55 (36) 



Ap r 



ruts 



V 



0.05 , (37) 



13 



4 Results 



4.1 A perfect detector 

We first find the sensitivity of the NC unitarity test for a perfect detector, i.e., a 
detector that can categorize each event correctly as CC muon or NC, with no mis- 
identification and 100% efficiency The figure of merit for a perfect detector is given 
by Eq. El with (ncnc = 1- We show the 3a sensitivity for P^ s (i.e., the minimum 
value of P^ s for which P^ s = 35P^ IS ) versus N® (the number of CC muons expected 
in the detector with no oscillations) for several values of the NC systematic error in 
Fig. [T] (the dotted curves). At low statistics the sensitivity is very poor, and for high 
statistics the sensitivity approaches the asymptotic limit of 3ejvc/(l + ^nc), where 
€nc = 3N NC /N NC is the fractional NC normalization uncertainty. 



4.2 More realistic K2K-like and T2K-like experiments 

Next we find the NC sensitivity for the K2K-like detector described in Sec. 13.11 for 
the case P^e ~ 0. We generated 400,000 neutrino events using the NEUGEN simu- 
lator, from which the normalized mis-identification factors f x>y were calculated. For 
a given set of probabilities P xy , the values of Rnc an d were calculated, and the 
corresponding measured value of P^ s was determined from Eq. ^2 The uncertainty 
on P^ s was calculated using Eq. HU assuming the uncertainties 5 Rnc an d SR^ are 
uncorrelated and add in quadrature. The 3a sensitivity for P MS is shown in Fig. ^ 
(solid curves) for various values of ejvc for the case P MS = 1 — P^ (all z/ M oscillating 
to u s ; we will consider cases with nonzero P^ T later). For both low and high statistics 
the K2K-like 3a sensitivity can be approximated by 

pmin _ 3(1 + f^Nc){&RNc/RNc) /„ R x 

f " 1 + 3(1 +f lhNC )(6R NC /RNc)' { ] 

which can be derived from Eq. EE where factors quadratic in the f x%y are ignored. 
Since / M) jvc — 0.08 for our K2K-like experiment, the NC sensitivity is at most about 
1.08 worse than that of the perfect detector for large numbers of events where the 
statistical uncertainty becomes negligible compared to the systematic uncertainty. At 
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low statistics the efficiency becomes important and the K2K-like performance will be 
more than 1.08 worse than a perfect detector. 

The K2K-like curves in Fig. [2 are plotted for the simple K2K signals in Table ^ 
The corresponding curves for the more aggressive K2K-like signals are very similar to 
the simple case; the improved efficiencies are partially compensated for by the slightly 
higher value of f^NC- Thus the result is fairly insensitive to the exact signal criteria 
used. 

We next consider the effects of nonzero P MT . If we assume P^ r = 1 — P MS (i.e., 
= Py,e = 0), the curves are very close to those of the perfect detector, since the 
dominant mis-identification term f^NC does not contribute to R^c when P m = 0. 
If both P w and P^ r are both nonzero (with P^ e ~ 0), the results will lie somewhere 
between the curves for K2K-like and the perfect detector. 

Finally, we consider nonzero P Me , in which case R e must also be measured and 
P^ s is determined using Eq. |23 As discussed in Sec. 12.3.11 if P^ e is of order 0.1 or 
less (as indicated by oscillation bounds such as from the CHOOZ reactor), and if 
the misidentification factors are also of order 0.1 or less, then this case reduces to 
that where the v e are ignored. We have verified this numerically for the K2K-like 
misidentification factors in Tabled 

In summary, the sensitivity of the K2K-like detector to the NC signal is only 
slightly worse than that of a perfect detector, with the dominant loss of sensitivity 
coming from the mis-identification of CC muon events as NC. For comparison, in 
Fig. Q we have also shown sensitivity curves for the T2K-like experiment with a 
Gaussian beam spread. Since f^NC — 0.25 in this case, the sensitivity is about a 
factor of 1.25/1.08 = 1.16 worse than for K2K. 

4.3 The MINOS-like detector 

For the MINOS-like case, we generated 320,000 neutrino events using the NEUGEN 
simulator, and calculated the corresponding mis- identification factors. The 3a sensi- 
tivity for P^ s was calculated as described above for the case P^ s = 1 — P MM ; the results 
are shown in Fig. ^ At low statistics, the MINOS-like experiment does better than 
the K2K-like and T2K-like experiments because of the higher NC efficiency, but at 
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high statistics it does worse because of the larger mis-identification factors. 
4.4 Exclusion limit when P^ s = 

If a 3cr signal for P^ s is not observed, then an exclusion limit (upper bound) for P^ s 
can then be obtained. The 90% CL. exclusion limit for P^ s is shown in Fig. |2] for a 
perfect detector (dotted curves), K2K-like with basic signals (solid curves), and T2K- 
like with basic signals (dashed curves). To model realistic oscillation probabilities we 
have assumed a three-neutrino model assuming the parameters Sm^ = 2.0 x 10~ 3 eV 2 , 
sin 2 2^ 2 3 = 1.0, and sin 2 2# 13 = 0.1. At high statistics the relative values of the 
exclusion limits are approximately proportional to (1 + f^^Nc)-, similar to the 3a 
sensitivity levels calculated previously. As was the case for the 3a sensitivity, the 
MINOS-like detector does better than the K2K-like and T2K-like detectors at low 
statistics, due to the higher NC efficiency, but not as well at high statistics due to 
larger misidentification factors. 

5 Summary 

At low statistics (< 1000 events), experiments with a larger NC efficiency, such as our 
MINOS-like example, tend to have better sensitivity to the sterile oscillation proba- 
bility P^ s . At high statistics, the sensitivity in the cases we considered is primarily 
limited by the systematic uncertainty in the NC rate, e^c, and the contamination of 
the NC signal from CC events, f^NC (and NC contamination from CC r events, 
fr,NC, above r threshold). The best anticipated €nc is of order a few per cent, so 
the best 3a sensitivity and 90% C.L. exclusion limits that can be expected for the 
sterile oscillation probability will be of order 0.10-0.15 (0.2-0.3 for the oscillation 
amplitude). The lowest contamination rates are realized for the K2K-like and T2K- 
like cases. Significant improvements in these sterile probability sensitivies or limits 
can only be achieved by lowering the uncertainty in NC cross sections or improving 
the event selection criteria, both of which could prove to be challenging but very 
worthwhile. 
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Figure 1: Assuming Gaussian statistical uncertainties, the 3a sensitivity for measuring 
P^ s versus N® for fixed values of the NC systematic error for a perfect detector 
(dotted curves, using Eq.ll7j). the K2K-like experiment with our basic signal definition 
(solid), the T2K-like experiment (dashed), and the MINOS-like experiment (dash- 
dotted). The number of NC events without oscillations is N^ c = 0.156iV° The 
systematic uncertainties 8N®/N® and 5N®/N® are assumed to be 2%, except when 
€nc = ^^nc/^nc = 0' i n w hich case they are 0. The arrows indicate the approximate 
statistical sensitivities expected for the K2K, T2K and MINOS experiments. 
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Figure 2: Assuming Gaussian statistical uncertainties, the 90% CL. exclusion limit 
for P^ s versus N® for fixed values of the NC systematic error for a perfect detector 
(dotted curves), the K2K-like experiment with basic signals (solid), the T2K-like 
experiment with basic signals (dashed), and the MINO-like experiment (dash-dotted). 
Other assumptions are the same as in Fig. ^ 
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